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The weak-boson fusion process is expected to provide crucial information on Higgs 
boson couplings at the Large Hadron Collider at CERN. The achievable statistical 
accuracy demands comparison with next-to-leading order QCD calculations, which are 
presented here in the form of a fully flexible parton Monte Carlo program. QCD 
corrections are determined for jet distributions and are shown to be modest, of order 
5 to 10% in most cases, but reaching 30% occasionally. Remaining scale uncertainties 
range from order 5% or less for distributions to below ±2% for the Higgs boson cross 
section in typical weak-boson fusion search regions. 



1 Introduction 



The weak-boson fusion (WBF) process, qQ^qQH, is expected to provide a copious 
sources of Higgs bosons in pp-collisions at the Large Hadron Collider (LHC) at CERN. It 
can be visualized (see Fig. |TJ) as the elastic scattering of two ( ant i-) quarks, mediated by 
t-channel W or Z-exchange, with the Higgs boson radiated off the weak-boson propagator. 
Together with gluon fusion, it represents the most promising production process for Higgs 
boson discovery [TJ |2j. Once the Higgs boson has been found and its mass determined, 
the measurement of its couplings to gauge bosons and fermions will be of main interest jS]. 
Here WBF will be of central importance since it allows for independent observation in the 
H — >tt [4 , H^WW jS], H— >77 [S| and H— > invisible J7| channels. This multitude of 
channels is crucial for separating the effects of different Higgs boson couplings. 

The WBF measurements can be performed at the LHC with statistical accuracies on 
cross sections times decay branching ratios, a ■ B, reaching 5 to 10% [3j. In order to extract 
Higgs boson coupling constants with this full statistical power, a theoretical prediction of 
the Standard Model (SM) production cross section with error well below 10% is required, 
and this clearly entails knowledge of the next-to-leading order (NLO) QCD corrections. 

For the total Higgs boson production cross section via WBF these NLO corrections 
have been available for a decade jH] and they are relatively small, with i^-factors around 
1.05 to 1.1. These modest if -factors are another reason for the importance of Higgs boson 
production via WBF: theoretical uncertainties will not limit the precision of the coupling 
measurements. This is in contrast to the dominant gluon fusion channel where the if-factor 
is larger than 2 and residual uncertainties of 10-20% remain, even after the 2-loop corrections 
have been evaluated P ITUj. 

In order to distinguish the WBF Higgs boson signal from backgrounds, stringent cuts are 
required on the Higgs boson decay products as well as on the two forward quark jets which 
are characteristic for WBF. Typical cuts have an acceptance of less than 25% of the starting 
value for a ■ B. The question then arises whether the if -factors and the scale dependence 
determined for the inclusive production cross section [8, are valid for the Higgs boson search 
region also. This is best addressed by implementing the one-loop QCD corrections in a fully 
flexible NLO parton-level Monte Carlo program. 

We are presently developing such programs for a collection of relevant WBF processes, 
of which Higgs boson production, in the narrow resonance approximation, is the simplest 
example. The purpose of this paper then is twofold. First we use the Higgs boson signal 
process as our example to discuss the generic features of NLO QCD corrections to WBF 
processes. We use the subtraction method of Catani and Seymour [TT] throughout. In 
Section |21 we describe the handling of real emission singularities. We give explicit formulas 
for the finite contributions which remain after factorization of the initial-state collinear 
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Figure 1: Feynman graphs contributing to qQ^qQH at (a) tree level and (b) 
including virtual corrections to the upper quark line. 

singularities and after cancellation of divergences produced by soft and collinear final-state 
gluons against the corresponding terms in the virtual corrections. 

This procedure yields a regularized Monte Carlo program which allows us to determine 
infrared safe observables at NLO. The main features of the program, numerical tests, and 
parameters to be used in the later phenomenological discussion are described in Sectional 
In Section |U we then use this tool to address our second objective, a discussion of the QCD 
radiative corrections as a function of jet observables. We determine the i^-factors and the 
residual scale uncertainties for distributions of the tagging jets which represent the scattered 
quarks in WBF. In addition, we quantify the cross section error induced by uncertainties in 
the determination of parton distribution functions (pdf's). Pdf errors and scale variations 
in the phase-space regions relevant for the Higgs boson search turn out to be quite small 
(approximately 4% when combined) and thus indicate the small theoretical uncertainties 
required for reliable coupling measurements. Conclusion are presented in Section 

2 Subtraction terms for soft and collinear radiation 

At lowest order, Higgs boson production via weak-boson fusion is represented by a single 
Feynman graph, like the one depicted in Fig. [TJa) for qQ^qQH. We use this particular pro- 
cess to describe the QCD radiative corrections. Generalization to crossed processes (q^q 
and/or Q-^Q) is straightforward. Strictly speaking, the single Feynman graph picture is 
valid for different quark flavors on the two fermion lines only. For identical flavors annihila- 
tion processes, like qq—>Z*^ZH with subsequent decay Z^-qq or similar WH production 
channels, contribute as well. For qq—>qqH or qq—xjqH the interchange of identical quarks 
in the initial or final state needs to be considered in principal. However, in the phase-space 
regions where WBF can be observed experimentally, with widely separated quark jets of very 
large invariant mass, the interference of these additional graphs is strongly suppressed by 
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Figure 2: Real emission contributions to Higgs boson production via weak-boson 
fusion. Corrections for the upper quark line only are shown: gluon radiation ((a) 
and (b)) and gluon initiated processes ((c) and (d)). 



large momentum transfer in the weak-boson propagators. Color suppression further makes 
these effects negligible. In the following we systematically neglect any identical fermion 
effects. 

At NLO, the vertex corrections of Fig. H^b) and the real emission diagrams of Fig. |2] 
must be included. Because of the color singlet nature of the exchanged weak boson, any 
interference terms between sub-amplitudes with gluons attached to both the upper and the 
lower quark lines vanish identically at order a s . Hence it is sufficient to consider radiative 
corrections to a single quark line only, which we here take as the upper one. Corrections 
to the lower fermion line are an exact copy. We denote the amplitude for the real emission 
process 

q(p a ) + Q(pb)-+g(pi) + q(p 2 ) + Q(ps) + H(P) (2.1) 

depicted in Fig. |2fa) and (b) as A4f = A4f(p a ,pi,P2] <?), where q = p% + P2 — p a is the four 
momentum of the virtual weak boson, V, of virtuality Q 2 = — q 2 . 

The 3-parton phase-space integral of |A^^| 2 suffers from soft and collinear divergences. 
They are absorbed in a single counter term, which, in the notation of Ref. ^Tj, contains the 
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where C F = I and M. q B = A^g(p a ,p2j ?) is the Born amplitude for the lowest order process 
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This choice continuously interpolates between the singularities due to final-state soft gluons 
(pi— »0 corresponding to x— >1 and z—*l), collinear final-state partons (pi||p2 resulting in 
Pi -p2-^0 or x—*l) and gluon emission collinear to the initial-state anti-quark (pi— >(1 — x)p a 
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and z— The subtracted real emission amplitude squared, |jM 
finite phase-space integral of the real parton emission cross section 

rl r-l 1 

a^ LO (qQ-^qQHg) = dx a dx b f q/p (x a , fi F ) f Q / p {x b ,ii F ) — d$ 4 (pi,p 2 ,P3, P]Pa + Pb) 
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where s = (p a + p b ) 2 is the center-of-mass energy. The functions Fj 6> and Fj 2 ^ define the 
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jet algorithm for 3-parton and 2-parton final states and we obviously need Fj —n-j 
singular limits discussed above, i.e. the jet algorithm (and all observables) must be infrared 
and collinear safe. Being finite, the phase-space integral of Eq. f!2.7|) is evaluated numerically 
in D = 4 dimensions. Similarly, for the gluon initiated process 



g{p a ) + Q{p b ) -> g{pi) + g(p 2 ) + Q{p 3 ) + H{P) , 

the singular behavior for g—>-qq splitting is absorbed into the singular counter term 
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where T F = | and }A B and M. q B denote the Born amplitudes for the leading-order (LO) 
processes q(p a ) + Q(pfc)-><?(p 2 ) + Qips) + H(P) and q(p a ) + Q{p b )^q{p2) + Q{Pz) + H(P), 
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respectively. The subtraction of l-M 9 ^ from the real emission amplitude squared leads to a 
contribution to the subtracted 3-parton cross section analogous to the one given in Eq. (|2.7p . 

The singular counter terms are integrated analytically, in D = 4 — 2e dimensions, over 
the phase space of the collinear and/or soft final-state parton. Integrating Eq. (|2.2|) yields 
the contribution (we are using the notation of Ref. [TT] ) 
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(2.10) 



We have regularized the divergences using dimensional reduction. If we had used conven- 
tional dimensional regularization we would have obtained a finite piece equal to (10 — 47r 2 /3). 
The 1/e 2 and 1/e divergences cancel against the poles of the virtual correction, depicted in 
Fig. IHb). For the case at hand, the virtual correction amplitude M.y is particularly simple, 
leading to the divergent interference term 



2Re[M v M* B ] = \M% 



5 \2 



2tt 
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Q' 2 



r(i + e) 



+ c 



virt 



(2.11) 



Here we have included the finite contribution of the virtual diagram which is proportional to 
the Born amplitude. In dimensional reduction this contribution is given by cvi r t = 7r 2 /3 — 7 
( c v i r t = 7r 2 /3 — 8 in conventional dimensional regularization). 

Summing together the contributions from Eq. ()2.10|) and Eq. (j2.11j) . we obtain the fol- 
lowing finite 2-parton contribution to the NLO cross section 

(T2 LO (qQ^qQH) = dx a dx b f q/p (x a , fi F ) f Q / p (x b , y. F ) — d$ 3 (p 2 ,p 3 , P)Va + Pb) 

J J ZiS 



(P2,P3) 



1 + 



2tt 



C f 9 - -TT 2 + c V i rt 



.(2.12) 



The two a s terms, at distinct renormalization scales fiR a and fipa,, correspond to virtual 
corrections to the upper and the lower fermion line in Fig. Q respectively, and we have 
anticipated the possibility of using different scales (like the virtuality of the attached weak 
boson V) for the QCD corrections to the two fermion lines. 

The remaining divergent piece of the integral of the counter terms in Eqs. (12. 2 j) and ()2.9|) 
is proportional to the P qq and P 9q splitting functions and disappears after renormalization 
of the parton distribution functions. The surviving finite collinear terms are given by 
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and similarly for quark initiated processes. Here the anti-quark function f?, p (x, [ip, [Ir) is 
given by 
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(2.14) 

(2.15) 
(2.16) 
(2.17) 
. (2.18) 



Note that c virt exactly cancels between the contributions from Eq. ()2.12|) and Eq. (|2.18j) . 
This fact will be used below to numerically test our program. 

The same kernels define the quark functions fq/ p (x, fi F , /ir), which appear with the Born 
amplitude M^ipa^V^ q) in the analog of Eq. ()2.13|) for the qQ^qQH processes. The gluon 
distribution f g / p (x,/if) thus appears twice, multiplying the Born amplitudes squared l-M^] 2 
and l-M^I 2 in the quark and anti-quark functions. These two terms correspond to the two 
terms in Eq. (|2.9|) . after the 1/e collinear divergences have been factorized into the NLO 
parton distributions. 

Formulae identical to the ones given above for corrections to the upper line in the dia- 
grams of Fig. El apply to the case where the gluon is attached to the lower line (with a <-» b, 
P2 Vz)- As for the renormalization scale /Jr in Eq. (|2.12j) . we distinguish between the two 
factorization scales that appear for the upper and lower quark lines, calling them \x Fa and 
/ipb, when needed. 

A second class of gluon initiated processes arises from crossing the final-state gluon and 
the initial-state quark Q in the Feynman graphs of Fig. Ufa) and (b). The resulting process 
can be described as gq—*qVH with the virtual weak boson V undergoing the hadronic decay 
V—>-QQ. Such contributions are part of the radiative corrections to qq-^VH, they are 
suppressed in the WBF search regions with their large dijet invariant mass, and we do not 
include them in our calculation. 
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3 The NLO parton Monte Carlo program 



The cross section contributions discussed above for the qQ^qQH process and cross- 
ing related channels have been implemented in a parton-level Monte Carlo program. The 
tree-level amplitudes are calculated numerically, using the helicity-amplitude formalism of 
Ref . [12] . The Monte Carlo integration is performed with a modified version of VEGAS jTS] • 

The subtraction method requires the evaluation of real-emission amplitudes and, simul- 
taneously, Born amplitudes at related phase-space points (see e.g. Eqs. (12. 7j) and (|2.13j) ). 
In order to speed up the program, the contributions from (J 3 LO and cr^raU are calculated 
in parallel, as part of the 3-parton phase-space integral. Since the phase-space element 
factorizes jTTj, 

f f 1 f 1 Q 2 

d$ 4 (p 1 ,p 2 ,P3,P;p a + Pb) = dx dzd$ 3 (p 2 ,p 3 ,P;xp a +p b ) , (3.1) 

J Jo Jo lbir z x 

we can rewrite the finite collinear term of Eq. (J2.13j) as 

a a i 

^coii (QQ^QQH) = dx a dx b — ■ ^d$ 4 (p 1 ,p 2 ,P3,P;Pa + Pb) 

Jo Jo 2(p a +p b y 

X ifg/ P ( X a, Vf)M X ) + fq/p( X a, Hf) \B(x) + C (x) 

D(xx a 



+ xfq/ p (xX a ,IJ F ) 



1 - XX a 



B x) 



xf Q / P (x b ,»F) \M%\ 2 F { ? (p 2 ,p 3 ) , (3.2) 

where x and z are determined as in Eqs. ()2.5|) and ([2.6)1 . Equation (|3.2j) allows for stringent 
consistency checks of our program, since we can determine the finite collinear cross section 
either as part of the 2-parton or as part of the 3-parton phase-space integral. For example, 
because of the cancellation of c virt mentioned below Eq. (|2.18j) . the final result cannot depend 
on its value. We have checked this independence numerically, at the 3 x 10~ 4 level. Another 
method to test the program is to determine the (anti-)quark functions fa/ p ( x , Hf, Hr) by 
numerical integration of Eq. ()2.14|) . to then compute the finite collinear cross section together 
with the Born cross section, and to compare with the results of Eq. ()3.2|) . For all phase-space 
regions considered, our numerical program passes this test, with relative deviations of less 
than 2 x 10 -4 of the total Higgs boson cross section, which is the level of the Monte Carlo 
error. 

As a final check we have compared our total Higgs boson cross section with previous 
analytical results [H] , as calculated with the program of Spira [H] . We find agreement at or 
below the 1 x 10~ 3 level which is inside the Monte Carlo accuracy for this comparison. 
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The cross sections to be presented below are based on CTEQ6M parton distributions ^2] 
with a s (Mz) = 0.118 for all NLO results and CTEQ6L1 parton distributions with a s (Mz) = 
0.130 for all leading order cross sections. For all Z-exchange contributions the 6-quark is 
included as an initial and/or final-state massless parton. The 6-quark contributions are 
quite small, however, affecting the Higgs boson production cross section at the 3% level 
only. We choose raz = 91.188 GeV, oiqed = 1/128.93 and the measured value of Gf as our 
electroweak input parameters from which we obtain raw = 79.96 GeV and sin 2 9w = 0.2310, 
using LO electroweak relations. In order to reconstruct jets from the final-state partons, the 
fcr-algorithm [TH] as described in Ref. [T7] is used, with resolution parameter D = 0.8. 



4 Tagging jet properties at NLO 

The defining feature of weak-boson fusion events at hadron colliders is the presence of 
two forward tagging jets, which, at LO, correspond to the two scattered quarks in the process 
qQ—>-qQH. Their observation, in addition to exploiting the properties of the Higgs boson 
decay products, is crucial for the suppression of backgrounds jU EJ El U\- The stringent 
acceptance requirements imply that tagging jet distributions must be known precisely for a 
reliable prediction of the SM Higgs signal rate. Comparison of the observed Higgs produc- 
tion rate with this SM cross section, within cuts, then allows us to determine Higgs boson 
couplings [3| and, thus, the theoretical error of the SM cross section directly feeds into the 
uncertainty of measured couplings. 

The NLO corrections to the Higgs boson cross section do not depend on the phase space 
of the Higgs boson decay products because the Higgs boson, as a scalar, does not induce 
any spin correlations. It is therefore sufficient to analyze tagging jet distributions to gain a 
reliable impression of the size and the uncertainties of higher order QCD corrections. Since 
search strategies depend on the decay mode considered and will evolve with time, we here 
consider generic weak-boson fusion cuts only. They are chosen, however, to give a good 
approximation of the cuts suggested for specific Higgs boson search channels at the LHC. 
The phase space dependence of the QCD corrections and uncertainties, within these cuts, 
should then provide a reasonably complete and reliable picture. 

Using the fc^-algorithm, we calculate the partonic cross sections for events with at least 
two hard jets, which are required to have 

p Ti >20GeV, |^| < 4.5. (4.1) 

Here yj denotes the rapidity of the (massive) jet momentum which is reconstructed as the 
four- vector sum of massless partons of pseudorapidity \r]\ < 5. The Higgs boson decay 
products (generically called "leptons" in the following) are required to fall between the two 
tagging jets in rapidity and they should be well observable. While an exact definition of 
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criteria for the Higgs boson decay products will depend on the channel considered, we here 
substitute such specific requirements by generating isotropic Higgs boson decay into two 
massless "leptons" (which represent t + t~ or 77 or bb final states) and require 

p Ti > 20 GeV , M < 2.5 , AR je > 0.6 , (4.2) 

where Rj£ denotes the jet-lepton separation in the rapidity-azimuthal angle plane. In addition 
the two "leptons" are required to fall between the two tagging jets in rapidity 

yj,min < Vh,2 ^ Uj,max ■ 

(4.3) 

We do not specifically require the two tagging jets to reside in opposite detector hemispheres 
for the present analysis. Note that no reduction due to branching ratios for specific final 
states is included in our calculation: the cross section without any cuts corresponds to the 
total Higgs boson production cross section by weak-boson fusion. 

At LO, the signal process has exactly two massless final-state quarks, which are identified 
as the tagging jets, provided they pass the /c T -algorithm and the cuts described above. At 
NLO these jets may be composed of two partons (recombination effect) or we may encounter 
three well-separated partons, which satisfy the cuts of Eq. (j4.1j) and would give rise to three- 
jet events. As with LHC data, a choice needs to be made for selecting the tagging jets in 
such a multijet situation. We consider here the following two possibilities: 

1) Define the tagging jets as the two highest pt jets in the event. This ensures that 
the tagging jets are part of the hard scattering event. We call this selection the u px- 
method" for choosing tagging jets. 

2) Define the tagging jets as the two highest energy jets in the event. This selection favors 
the very energetic forward jets which are typical for weak-boson fusion processes. We 
call this selection the ".E-method" for choosing tagging jets. 

Backgrounds to weak-boson fusion are significantly suppressed by requiring a large rapidity 
separation of the two tagging jets. As a final cut, we require 

&U3j = \Vh-Vh\ > 4 , ( 4 - 4 ) 
which will be called the "rapidity gap cut" in the following. 

Cross sections, within the cuts of Eqs. (|4.1|) - (|4.4j) . are shown in Fig. OK a), as a function 
of the Higgs boson mass, m#. As for the total WBF cross section, the NLO effects are 
modest for the cross section within cuts, amounting to a 3-5% increase for the py-method 
of selecting tagging jets (solid red) and a 6-9% increase when the E'-method is used 1 . These 

1 The larger cross section for the _E-mcthod is due to events with a fairly energetic extra central jet. A 
veto on central jets of pxj > 20 GeV and rapidity between the two tagging jets, as suggested for the WBF 
selection, lowers the NLO cross section to 0.97 x olo for the pr-method and 0.93 x olo for the i?-method. 
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Figure 3: Effect of QCD radiative corrections on the Higgs boson production cross 
section via WBF, as a function of the Higgs boson mass, run- Results are given 
at LO (black dotted) and at NLO for the pT-method (solid red) and the E-method 
(dashed blue) for defining tagging jets. Panel (a) gives the total cross section within 
the cuts of Eqs. \4-l\) - \4-4\l - The corresponding scale dependence, for variation of 
fjLji and jip by a factor of 2, is shown in panel (b). See text for details. 



.fT-factors, and their scale dependence, are shown in Fig. E^b). Here the i^-factor is defined 

as 

<t™(ji f = q.) ' ^> 
i.e. the cross section is normalized to the LO cross section, determined with CTEQ6L1 parton 
distributions, and a factorization scale which is set to the virtuality of the weak boson which 
is attached to a given quark line. 

We have investigated two general scale choices. First we consider the Higgs boson mass 
as the relevant hard scale, i.e. we set 

Hf = i F m H , fji R = £ R m H . (4.6) 

As a second option we consider the virtuality of the exchanged weak boson. Specifically, 
independent scales Qi are determined for radiative correction on the upper and the lower 
quark line, and we set 

fJ-Fi = &Qi , = £,nQi ■ (4.7) 
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Figure 4: Variation of the total cross section, within cuts, due to errors in the 
parton distribution functions, as a function of mn- The central solid line corre- 
sponds to the "best fit" CTEQ6M pdf, while the upper and lower curves define the 
pdf error band, which is determined from the 40 error eigenvectors in the CTEQ6M 
set (CTEQ6M101-CTEQ6M140), adding cross section deviations in quadrature. 

This choice is motivated by the picture of WBF as two independent DIS events, with inde- 
pendent radiative corrections on the two electroweak boson vertices. In general we find the 
largest scale variations when we vary the renormalization scale and the factorization scale in 
the same direction. We only show results for this case, £ = = £p, in the following. The 
curves in Fig. Efb) correspond to the largest variations found for £ = 1/2 and £ = 2 when 
considering both scale choices simultaneously. The residual scale uncertainty is about ±5% 
at LO and reduces to below ±2% at NLO. 

In addition to missing higher order corrections, the theoretical error of the WBF cross 
section is dominated by uncertainties in the determination of the parton distribution func- 
tions. We have investigated this dependence by calculating the total Higgs boson cross 
section, within the cuts of Eqs. (|Q ]l -(fO Jl . for the 40 pdf's in the CTEQ6Mxxx (xxx = 
101-140) set. They correspond to extremal plus/minus variations in the directions of the 
20 error eigenvectors of the Hessian of the CTEQ6M fitting parameters [15] . Adding the 
maximum deviations for each error eigenvector in quadrature, one obtains the blue dashed 
lines in Fig. which define the pdf error band. We find a uniform ±3.5% pdf uncertainty 
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Figure 5: Transverse momentum distribution of the softer tagging jet for the the 
PT-method (solid red) and the E-method (dashed blue) of defining tagging jets, for 
mu = 120 GeV. The right-hand panels give the K -factors (black dash-dotted line) 
and the scale variation of the NLO results. Solid colored curves correspond to [if = 
1*>R = £,Qi an d dashed colored curves are for fip = /j>r = with £ = 1/2 and 2. 
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1.3 




Figure 6: Higgs boson production cross section as a function of the smaller of 
the absolute value of the two tagging jet rapidities, da/d\y\™g (in fb, for ran = 
120 GeV). Results are shown at LO (dotted black) and at NLO for the pT-method 
(solid red) and the E-method (dashed blue) of defining tagging jets. The right-hand 
panel gives the K -factor (black dash-dotted line) and the scale variation of the NLO 
result for the E-method. Colored curves for the scale dependence are labeled as in 
Fig.\% 



of the total cross section over the entire range of m# shown. 

Scale and pdf uncertainties exhibit little dependence on the Higgs boson mass. We there- 
fore limit our investigation to a single, representative Higgs boson mass for the remaining 
discussion, which we take as m# = 120 GeV. 

While the scale dependence of the integrated Higgs boson production cross section is 
quite weak, the same need not be true for the shape of distributions which will be used to 
discriminate between Higgs boson signal and various backgrounds. Having a fully flexible 
NLO Monte Carlo program at hand, we can investigate this question. Crucial distributions 
for the detection efficiency of the signal are the transverse momentum and the rapidity of 
the tagging jets. In Fig. El the cross section is shown as a function of p™™ ag , the smaller of 
the two tagging jet transverse momenta. At LO, the tagging jets are uniquely defined, but 
at NLO one finds relatively large differences between the pr-method (solid red curves in the 
top panels) and the .E-method (dashed blue curves in bottom panels). The right-hand-side 
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Figure 7: Rapidity separation of the two tagging jets for mn = 120 GeV. In the 
left-hand panel, da /d Ayjj (in fb) is shown at LO (dotted black) and at NLO (solid 
red), for the px-method of defining tagging jets. The right-hand panel gives the 
corresponding K-factor (black dash-dotted line) and the scale variation of the NLO 
results. Colored curves for the scale dependence are labeled as in Fig. 



panels give the corresponding .fT-factors, as defined in Eq. (J4.5j) . (black dash-dotted lines) 
and the ratio of NLO differential distributions for different scale choices. Shown are the 
ratios 

da NLO { f M F = f iR = Q i ) 

for \i = 2 ±1 Q i (solid lines) and fi = 2 ±1 m H (dashed lines). While the .fT-factor is modest 
for the pt~ method, it reaches values around 1.3 in the threshold region for the .E-method. 
This strong rise at NLO is due to hard forward gluon jets being misidentified as tagging jets 
in the E'-method. This problem was recognized previously in parton shower Monte Carlo 
simulations and has prompted a preference for the pr-method ^H]- I n spite of the large 
ET-factor, however, the residual scale uncertainty is small, ranging from -4% to +2% for the 
Pr-method and -2% to +5% for the E-method. 

The more forward selection of tagging jets in the E-method is most obvious in the rapidity 
distributions of Figs. E]and[7| In Fig. El the rapidity of the more central of the tagging jets, 
|y|™ g n , is shown. At NLO, the tagging jets are slightly more forward than at tree level, leading 
to a i^-factor which varies appreciably over phase space. This |y|™ g n - dependence is shown in 
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Figure 8: Transverse momentum distribution of the harder tagging jet, for mu = 
120 GeV. In the left-hand panel, da/dp^^ (in fb/GeV) is shown at LO (dotted 
black) and at NLO (solid red), for the px-method of defining tagging jets. The right- 
hand panel gives the corresponding K -factor (black dash-dotted line) and the scale 
variation of the NLO results. Colored curves for the scale dependence are labeled as 
in Fig. El 



the right-hand panel for the .E-method, together with the residual scale dependence at NLO. 
Again, scale variations of less than ±4% are found over virtually the entire phase space. For 
the pr-method, similar scale variations arise, as shown in Fig. [7|for the rapidity separation 
between the two tagging jets, where the cuts of Eqs. (j4.1|) - ()4.3|) have been imposed. Figured 
demonstrates that the wide separation of the tagging jets, which is important for rejection 
of QCD backgrounds, does survive at NLO. In fact, the tagging jet separation even increases 
slightly, making a separation cut like Ay = \yj 1 — yj 2 \ > 4 even more effective than at LO. 

In all distributions considered so far, no clear preference emerges on whether to choose the 
weak-boson virtuality, Qi, or m H as the hard scale. While both choices are acceptable, the 
transverse momentum distributions show somewhat smaller scale variations for [i = £,Qi than 
H = £m#. The effect is most pronounced in the high px tail of the tagging jet distributions. 
When considering da/dp™^ as shown in Fig. |HJ the scale variation increases to +10% at 
large pt when /i = is taken, while the same distribution for fi = stays in a narrow 
±2% band. This observation provides another reason for our default scale choice, /i = Qi. 
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Unlike the tagging jets considered so far, distributions of the Higgs decay products show 
little change in shape at NLO. 

5 Conclusions 

Weak-boson fusion processes will play an important role at future hadron colliders, most 
notably as a probe for electroweak symmetry breaking. For the particular case of Higgs boson 
production, we have presented a first analysis of the size and of the remaining uncertainties 
of NLO QCD corrections to jet distributions in WBF. 

As for the inclusive WBF cross section, QCD corrections to distributions are of modest 
size, of order 10%, but occasionally they reach larger values. These corrections are strongly 
phase-space dependent for jet observables and an overall if-factor multiplying the LO dis- 
tributions is not an adequate approximation. Within the phase-space region relevant for 
Higgs boson searches, we find differential i^-factors as small as 0.9 or as large as 1.3. These 
corrections need to be taken into account for Higgs coupling measurements, and our NLO 
Monte Carlo program, or the recently released analogous program in the MCFM package JH] , 
provide the necessary tools. 

After inclusion of the one-loop QCD corrections, remaining uncertainties due to as yet 
uncalculated higher order terms, can be estimated by considering scale variations of the 
NLO cross section. Using the Higgs boson mass, m H , and the weak-boson virtuality, Q iy 
as potential hard scales, we find that these remaining scale dependencies are quite small. 
Varying renormalization or factorization scales by a factor of two away from these two central 
values results in typical changes of the NLO differential cross sections by ±2% or less. The 
uncertainty bands for \i = ^m H and fi = ^Q, t typically overlap, yielding combined scale 
uncertainties of less than ±3% in most cases, occasionally rising up to order 5% at the 
edges of phase space. Moreover, the variation in different regions typically cancels in the 
integrated Higgs cross section, within cuts, leading to uncertainties due to higher order effects 
of ±2% (see Fig. |3J), even when considering different hard scales. The remaining theoretical 
uncertainty on the measurable Higgs cross section, thus, is well below expected statistical 
errors, except for the H^WW search for Higgs masses around 170 GeV, where the high 
LHC rate allows statistical errors as low as 3%. In addition, pdf uncertainties for the total 
cross section are of order ±3.5% over the range 100 GeV < mu < 200 GeV. This means 
that the SM Higgs boson production cross section via WBF can be predicted, at present, 
with a theoretical error of about ±4%. 

The expected size of the LHC Higgs signal is enhanced slightly by the NLO QCD correc- 
tions. In addition to a i^-factor slightly above unity this is due to a small shift of the tagging 
jets to higher rapidities, still well inside the detector coverage but moving the tagging jets 



17 



slightly farther apart and hence allowing a better differentiation from QCD backgrounds. 

The techniques described in this paper work in a very similar fashion for other weak- 
boson fusion processes. We are planning to extend our work to include also W and Z boson 
production via weak-boson fusion, and to make these programs generally available. 
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